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Abstract 

We demonstrate that by carefuhy analyzing the temperature dependent characteristics of the 
FV measurements for a given complex system it is possible to determine whether it is composed of 
a single, double or multiple quantum-dot structure. Our approach is based on the orthodox theory 
for a double-dot case and is capable of simulating I-V characteristics of systems with any resistance 
and capacitance values and for temperatures corresponding to thermal energies larger than the dot 
level spacing. We compare FV characteristics of single-dot and double-dot systems and show that 
for a given measured FV curve considering the possibility of a second dot is equivalent to decreasing 
the fit temperature. Thus, our method allows one to gain information about the structure of an 
experimental system based on an FV measurement. 

PACS numbers: 73.63.Kv; 73.23.Hk 
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Much of the study of charging effects in quantum systems has focused on a single-dot 
system, in which a metaUic island is coupled to two metallic leads via tunnel junctions 
If the electron density is large enough so that the discreet energy level spacing is 
negligible compared to the other energy scales, the relevant scales are the thermal energy, 
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Ex = i^sT, and the charging energy, Ec = |t, where C is the capacitance of the dot. For 
low temperatures and small bias {Et < eV < Ec), no current can flow through the dot 
leading to a Coulomb blockade in the I-V characteristic. An effective method to treat such a 
system is the orthodox theory . In this framework the quantum dot is represented by a 
Double Barrier Tunneling Junction (DBTJ) in which each junction is modeled as a resistor 
and a capacitor connected in parallel. For a given voltage, V, a distribution function, p, 
determines the probability at time t for the system to have N extra electrons on the island. 
Solving the master equation for p enables one to derive the current- volt age relation as a 
function of the tunneling rates. The orthodox model assumes that the tunneling events are 
sequential and an electron loses its phase coherence during a tunnel process, thus, quantum 
corrections are not taken into account P,^]- 

This approach has been very successful in analyzing the behavior of a single quantum 
dot. In many experimental systems, however, the exact structure is unknown. Even if a 
sample exhibits Coulomb-blockade-like features, one can not always be certain that only 
one simple quantum dot is involved in the transport process. An example for such a system 
is demonstrated in figure 1 which shows the I-V characteristic of a Co nanowire, 10pm 
long and 200nm wide j^. The wire was evaporated on a step-edge structure and was 
allowed to oxidize in atmosphere. It is seen that the I-V curve shows Coulomb-blockade- 
like behavior, presumably, due to the oxidation that gives rise to the formation of metallic 
islands separated by nanoconstrictions. However, a-priori it is impossible to know whether 
the structure consists of one dominant quantum dot, two dots or perhaps even more. A 
similar situation occurs in many experimental configurations and often it is desirable to find 
a way to determine the exact structure of the sample. In this paper we show that it is, 
in principle, possible to fit given I-V curves to both single and double-dot configurations 
using the orthodox model when the capacitance and resistance of the barriers as well as the 
temperature are treated as fit parameters. Nevertheless, since the temperature is usually 
well controlled in the experimental set up, it is possible to determine whether a system 
contains one or more quantum dots based on I-V measurements at a given temperature. 
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FIG. 1: I-V characteristic at T=10K of a Co wire evaporated on a step-edge. The inset shows an 
SEM picture of the wire (white hne) grown between two large Co electrodes (gray pads). 



Most research studying the transport through double-dot system have taken into account 
the discrete energy spectrum in the dots 0,17,0]. In our approach, the more common case 
for metallic dots, i.e. discrete energy level spacing much smaller than KbT, is considered. 

The orthodox theory solution to the single dot (DBTJ) case is well known . Here 
we shall describe the double-dot case which contains three tunnel junctions and is named 
Triple Barrier Tunnel Junction (TBTJ). As in the DBTJ case, we assume that the tunnel- 
ing coefficient is low and tunneling events are sequential and non-coherent, hence quantum 
interference effects are neglected. In our TBTJ model each of the three junctions i=l,2,3 is 
characterized by a tunneling resistance Ri and a capacitance Ci (see Fig. 2). The "state" of 
the system is determined by the number of excess electrons on each grain (^1,722). Similar 
to the case of DBTJ for T=0 the distribution function is sharply peaked for each voltage 
around a most probable state {nl,n2)- However, in contrast to the DBTJ case, one cannot 
determine this state analytically and a more com plex method is required. Previous calcula- 
tions on I-V characteristics of TBTJ systems P. Iiol] overcame this difficulty by restricting 
the junction parameters (the resistance between the dots was assumed to be much higher 
dot-lead resistances), and the temperature was taken to be zero. We suggest an approach 
that allows one to calculate the distribution function in a general case without making any 
assumptions on the system. Moreover we calculate the probability value for any state, thus 
we can simulate the I-V curves for any given temperature. 
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FIG. 2: A schematic illustration of a system through which three electrons are allowed to pass. 
Each box is divided to two parts where each denotes one dot. Each dot is allowed to contain 0-3 
electrons. The arrows indicate a permitted tunneling between states. For the current calculation 
one has to sum over the transitions for a specific kind of arrow which are related to a specific 
junction: middle barrier (solid line), left barrier (dashed line) and right barrier (dotted line). The 
insert shows a schematic representation of a system composed of three tunnel junctions coupled in 
series. Each junction is characterized by a capacitance and resistance. 

Applying voltage to a TBTJ system causes three different voltage drops Vi across the 
tunnel junctions i: 



^ eC2C3V 

{C1C2 + C2C3 + C3C1) 

V, = '^'^'^ (1) 

(C*lC*2 + C2C3 + C3C1) 

^ eC,C2V 

{C1C2 + C2C3 + C3C1) 

Accordingly, six tunneling rates have to be considered. Each tunneling rate depends 
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on the energy difference before and after the tunnehng event and on the resistance of the 
relevant junction. The tunnehng can be derived from Fermi's Golden rule: 



±^ e^R,(l-exp{-A{El,)/k^T)y 
where T^^^ is the electron tunnehng rate on (+) or off (-) dot k across junction i. The energy 
differences, E!^^ can be derived by subtracting the electrostatic energy on the island after 
the tunneling event from that before the tunneling adding the gain in energy due to the 
voltage drop. The total energy differences due to tunneling of electrons from the dots to the 
leads are given by: 



.pi _ {C2 + C,)[iQ^±ef-Qj] 

- 2{C,C, + C,C, + CO ^ 

where Q\ — erii, Q2 — en2- 

Taking into account the electrostatic energy differences of both the dots the energy dif- 
ference due to tunneling of electrons from one dot to another is given by: 



2(CiC2 + C2C3 + C3C1] 



^ {C, + C2)[{Q2±ey-Ql] ^ 
2(CiC2 + C2C3 + C3C1) 



(5) 



For the TBTJ p is the probability to have ni and ^2 electrons on the first and second 
grain respectively, hence the master equation in a double-dot system is: 

dp{V,ni,n2,t) , . 

— = p{v, Til - 1, n2, t)r+^ + p{v, Til + 1, n2, t)r_^ (6) 

+p{v, m, 712 - 1, t)r% + p{v, ni, 712 + 1, ^)ri2 
+p{v, m - 1, 7i2 + 1, t)rl^ + p{v, 7ii + 1, 7i2 - 1, t)r\ 

-p{V, n,,n2, t)[Tl, + T'_, + F^^ + F^ + F^^ + ^-2]- 

In order to find the steady state solution of the distribution function we take ^ = 0. 
Solving this equation requires a constraint on the number of electrons permitted to pass 



through the system. For a specific number of electrons rie, there exists a specific number 
of states, N. In Fig. 2 we show a schematic drawing of the transitions between the charge 
configuration of a double-dot system allowing the addition of up to three electrons to the 
system. 

For each state we manipulate the master equation, thus achieving a system of N linear 
equations, where N is the number of states. This system is described by the formula: 

r-p = o, (7) 

where F is a rate matrix and p is the states vector. The sum over all the probabilities should 
be one. For simplicity we add the normalization condition in the last row of the rate matrix. 
Thus, the previous equation takes the form: 



-(^1 + ^2) 



^ -2 



^ +2 

^ +1 



-(n2 + ni + ni + ri,) 



+2) 





< 0,0 > 






< 0,1 > 






< 1,0 > 






< ne, > 





(8) 

By solving this numerically we obtain the distribution function and the vector p. The 
current for a given voltage is derived by summing over all the possibilities for an electron to 
pass through a certain junction: 



I{V)= J2 C±,-p{V,n,,n2)-Tl,{n,,n2), (9) 

k,i,ni,n2 

where C^^, = ± is determined by the direction of the tunneling. 

The results presented in this paper are obtained for a system of up to 4 electrons tunneling 
through the dot. In order to verify that this docs not lead to a considerable loss of information 
we plot the numerical I-V characteristics of a typical TBTJ for 1 to 4 tunneling electrons 
in Fig. 3. It can be seen that the curves for 3 and 4 electrons practically coincide for the 
relevant voltage regime, hence we conclude that further increasing the number of electrons 
would not have a noticeable effect on the I-V characteristics. 
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FIG. 3: I-V curves for the TBTJ model T=OK, Ci = C2 = C3 = 5 • lO'^^F, = R2 = R3 = 
1 • 10^0. Each number in the graph indicates the number of electrons passing through the system. 
Curves 3 and 4 which corresponds to three and four electrons are, for all practical purposes, 
identical (at least in the range of the Coulomb Blockade). 

Once more than two barriers are considered the system's parameter space (i.e., possible 
values of Ri, Ci) becomes large. It is helpful then to gain some information out of the 
general properties of the I-V curve. In this paper we shall concentrate on the simplest case 
of symmetric I-V curves with no prominent staircase features, which are surprisingly common 
in the experiments that are discussed later. Since in the Coulomb blockade range the case 
for which Ri, Ci differs much from R3, C3 results in a non-symmetric I-V curve and in some 
cases leads to staircases with different widths and heights, we will only consider cases for 
which Ri,Ci is similar to R^^C^. Another crucial factor for the I-V characteristics are the 
parameters of the middle barrier {R2,C2)- Choosing the parameters of the middle barrier 
significantly different than the other barriers results in I-V curves that show more prominent 
staircase structure. In Fig. 4. we compare the I-V curves for the three cases discussed above. 
The first case in which the parameters of the barriers are significantly different shows a 
pronounced non-symmetric I-V curve. In the second case, for which Ri,Ci = R^.C^ ^ 
R21C21 the I-V curve shows prominent staircase jumps. On the other hand for the case in 
which all barriers are equal the I-V curve is symmetric and smooth. 
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FIG. 4: Theoretical I-V curves at T=Ok for three cases. Curve 1 is the curve obtained for a pure 
symmetric system with the parameters Ci = C2 = C3 = 5 • W~^^F, Ri = R2 = R3 = 1 ■ lO^O. 
Curve 2 was obtained on a system in which the parameters of middle barrier differ from rest 
barriers: Ci = C3 = 7 • IQ-'^^F, Ri = R^ = 1 ■ 10^ C2 = 3 • 10-l^ i?2 = 4 • lO^O. This results in 
prominent staircases. Curve 3 shows non-symmetric behavior where the dot-lead barriers are not 
equal: Ci = C2 = 5 • IQ-^^F, R^ = R^ = I ■ 10^, C3 = 2 ■ lO-^^O, Rz = 2- lO^O 



We demonstrate the effectiveness of our analysis in determining the structure of a complex 
system by applying it to Co nanowires such as that depicted in Fig. 1. Fig. 5 shows the 
experimental results and the numerical fits of I-V characteristics of a typical nanowire taken 
at different temperatures. We used the orthodox theory to fit these data using DBTJ and 
TBTJ models. For the DBTJ we were able to obtain reasonable fits only using much higher 
temperatures than those of the experiment. Moreover, the ratio between the measured and 
calculated temperatures increased with T. Using our TBTJ calculation we were able to fit 



the data using the correct measurement temperatures jl^]. This remarkable agreement for 
the different temperatures reinforces our confidence in the validity of the two-dot model to 
this experimental system. The diameters of the metallic islands according to our fit are 
found to be ~ 30nm. This is a reasonable value since it is close to the width of the wire. 

We have applied this analysis to other wires. In some cases, even the TBTJ model yields 
good fits only for temperatures much higher than the experimental T. We suspect that these 
samples contain a larger number of islands for which a model that takes into account four 
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FIG. 5: DBTJ (right) and TBTJ (left) theoretical fits (dashed lines) to experimental I-V experi- 
mental curves (solid lines) for Co nanowire samples measured at lOK (A), 20K (B), and 30K (C). 
The fitting parameters for DBTJ: = R2 = 7 ■ lO^f), Ci = C2 = 4 • lO'^^F, T=50K (A), 65K 
(B), 90K (C), for TBTJ: = R2 = = A ■ lO^J], Ci = C2 = C-^ = 1.55 • IQ-^^F, T=10K (A), 
20K (B), 30K (C). 

or more tunnel junctions is required. 

Another complex system in which this analysis method may prove useful is a disordered 
granular sample. In such a system the geometrical structure may be known but the precise 
elements that dominate the transport are unidentified. We applied our analysis to a 400nm 
sample of 20— 40nm grains placed on an insulating substrate and separated by a few nm of an 
insulating matrix. Though the sample contains about 400 grains, not all of them participate 
in the transport due to the hopping nature of the electric conductivity. Actually, it has been 



shown llf that for low temperatures the transport is dominated by hopping through one or 
two grains. I-V characteristics of such a system show highly non-ohmic behavior that can 
be interpreted as signs for charging effects. Fig. 6. shows an I-V curve taken at T = 
together with fits to DBTJ and TBTJ models. Again, we were not able to fit the results 
to a single dot system without increasing the fitting temperature considerably. Using the 
TBTJ model, on the other hand, we were able to obtain a very good fit for T = QK. From 
the fitting parameters we find that the diameter of the grains is about 40 nm, which is in 
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FIG. 6: Experimental results (solid lines) and theoretical fits (dashed lines) for the I-V characteris- 
tics of the granular systems. The left panel is the fit for the DBTJ at T = 30K,Ci = C2 = 
R^ = R^ = 1.55 • lO^n. The right panel is for the TBTJ at T = 6K, Ci = C2 = C3 = 2.8 ■ IQ-^^F, 
R^ = R^ = R^ = l. lO^Q 

accordance with AFM pictures obtained on these systems. 

Both examples demonstrate that for a given experimental I-V curve, using the TBTJ 
model had a similar effect to that of using a DBTJ with higher temperature. This can be 
expected since the voltage drop in the Coulomb blockade regime is divided to the contribu- 
tions of the two dots, each one contributing an energy which has to be compared to KbT. 
Thus, the measurement temperature can be an important tool for determining the number 
of dots in a complex Coulomb blockade system. For further confirmation of the analysis it 
is recommendable to acquire I-V curves for different temperatures and repeat the procedure 
as demonstrated in Fig. 5. 

As previously discussed, since the typical experimental I-V curves considered in this pa- 
per were symmetric and had no pronounced staircases, we choose all barrier parameters to 
be equal in the fitting procedure. The assumption that all i?jS and CjS in the experimental 
structures are identical is clearly unrealistic. Nevertheless, we find that a finite tempera- 
ture smears the difference between barriers which exhibit similar parameters. The effect of 
temperature is demonstrated in Fig. 7 where the TBTJ I-V curves for two different values 
of the middle barrier parameters (i?2, C2 = Ri, Ci = -R3, C3 and R2, C2 ^ Ri,Ci = R3, C3) 
at zero temperature and at T = 6K are presented. At zero temperature there is a clear 
staircase structure for the latter case. At T = 6K, on the other hand, the staircase struc- 
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FIG. 7: Theoretical results for TBTJ at T=OK (solid lines) and at T=6K (dotted lines). The left 
panel is at Ci = C3 = 2.8 • 10'^^ F, i?i = i?3 = 5 • W^Q, R2 = 2 ■ IQ^Q, C2 = 2 ■ IO-17. The right 
panel is for the total symmetric case: Ci = C2 = C3 = 2.8 • 10~^^F, Ri = R2 = R3 = 1 ■ lO^O. 
One can notice that the staircases in the left panel are more prominent than in right panel and the 
dotted lines at both cases smears out the staircases.) 

ture is smeared and the I-V curve is similar for both identical and non-identical barriers. 
Hence, measured I-V curves cannot determine the precise parameters of the barriers for the 
experimental realizations depicted in Figs. 5 and 6 and the best we can do is to estimate 
the barrier parameters to be roughly equal. For a more exact evaluation of the barrier 
parameters additional measurements at lower temperatures are required. 

In summary, we propose that by analyzing the I-V characteristics one may identify sys- 
tems which are more complicated than the conventional single-dot-double-barrier system 
although no former knowledge of the system is assumed. This enables to determine the 
geometrical structure of a complex quantum system. In the current work we have imple- 
mented our approach to identifying two-dot-triple-barrier configuration. Future work should 
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focus on extending this method to apply to more complex systems and, eventually, to find a 
way to determine the precise number of dots in an experimental system. This research was 
supported by the Israeli Science Foundation (grants number 326/02-3, 877/04). 
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